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We derive compact analytical expressions for the coupled spin-charge susceptibility of a clean 
helical metal at the surface of a three dimensional topological insulator (TI). These expressions lead 
to unconventional non-coUinear RKKY interactions between two impurity magnetic moments placed 
on the surface of a TI, and predict the generation of electric currents by time-dependent magnetic 
moments. We determine the influence of gate and bias voltages on the interlayer exchange coupling 
between two single-domain ferromagnetic monolayers deposited on top of a TI. 
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Introduction. — The discovery of three dimensional 
topological insulators^ has unveiled a promising avenue 
for low-dissipation spintronics applications. The sur- 
faces of these materials host topologically protected, gap- 
less, strongly spin-orbit (SO) coupled conducting states 
dubbed helical metals. These are endowed with novel 
magnetoelectric properties when time-reversal (7~) sym- 
metry is broken^Ti^. Albeit less exotic, the magnetoelec- 
tric response of 7"-invariant helical metals is also pecu- 
liar and potentially useful. Recent theoretical studies^"— 
have discussed the spin-spin reponse of T-invariant he- 
lical metals, with an emphasis on RKKY interactions. 
This paper complements and extends previous work by 
presenting fully analytical expressions for the spin-spin as 
well as spin-charge response of a clean surface. In antici- 
pation of future experiments we predict that (i) the heli- 
cal metal mediates Dzyaloshinskii-Moriya coupling when 
the chemical potential lies away from the neutral point, 
(ii) electric fields may be used to generate and tune un- 
conventional RKKY interactions; (iii) precessing mag- 
netic moments produce alternating electric currents. 

Spin-charge response function. — The low-energy effec- 
tive Hamiltonian describing a helical metal living on the 
surface of a three dimensional TI isi 



"Ho = —iv 



a/3 



where v is the Fermi velocity, cr = {a^ , ,a^) are Pauli 
matrices denoting real spin, a,/3 G {t:4-}i r = (s;,?/), z 
is a unit vector normal to the surface and h = 1. The 
density matrix encoding the ground state charge and spin 
densities is 
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of frequency and wave vector q — q(cos0, sine 
density matrix changes via 
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where 
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is the dynamical spin-charge susceptibility and 
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is the Green's function for the surface statesi^i is the 
step function and /i ~ vkp is the chemical potential. 
From Eqs. (O and (O we arrive at 
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the form of which is compatible with symmetry ar- 
guments discussed in the context of ordinary two- 
dimensional electron systems (2DES) with Rashba SO 
interactions^. The charge-charge response function xoo 
is identical to that of grapheneiS and there is no cou- 
pUng between the charge density and the ^-component 
of the spin density. Adopting the formalism explained 
in Ref. [ll| we have derived the following concise expres- 
sions for the coefficients {/i}: 



where a G {0, x, y, z} labels charge (0) or spin (a;, y, z) 
sectors and cr*' = 1 is the identity matrix. Under a weak 
external potential 
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with X = 2kp/q and 
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with y = {2kF + fl/v)/q. We omitted 7^ expressions 
for /i,2.3; they are cumbersome and will not be needed 
below. /i(0, 0) = /2(0, 0) = indicates that, unlike in 
ordinary 2DES with SO interactions^^, a uniform and 
static in-plane magnetic field does not spin-polarize the 
helical metal: it merely shifts the location of the Dirac 
cone in momentum space^^. In contrast, uniform and 
static perpendicular magnetic fields elicit a paramagnetic 
response oc /3(0, 0) 0. Finally, /o(q, 0) = means that 
charge and spin sectors decouple in the static limit. 

RKKY interaction in electric equilibrium — Consider 
two localized spins Si and S2 placed at Ri and R2 on 
the surface of a 3D TI. Their indirect RKKY— coupling 
mediated by the helical metal is 



RKKY 
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where i,j G {x,y,z}, Jij = -A^jXij(R-): R = Ri - R2 
and Xij(R) is the spatial Fourier transform of Xij(qiO). 
In addition Xij are the exchange integrals (in units of 
energy x area) between the magnetic impurities and the 
surface states. Invoking spin rotational invariance in the 
xy plane we have Xj,^ = Xyy = X^y = Xy^ ^ X\\, X^z ^ X± 

and Xxz = Xzx = Xyz = X^y = ^ A|| Aj_. 
Combining Eqs. ([7]) and ([8]) we arrive at 



i^RKKY = A{SlSl^S\Sl) 
+ C(Si •R)(S2 -R) 



BSISI (11) 
^i^[Rx (Si XS2)]. 



which agrees with the results derived independently in 
Ref. Q . We are able to extract analytical expressions for 
the RKKY coefficients in limiting cases: 
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for kpR » 1 (see Fig[Bt)- In Eq. (fT2l) we have followed 
the prescription of Ref. [l5| to regularize ultraviolet diver- 
gences. In the derivation of Eq. we have integrated 



by parts and have used asymptotic expansions for inte- 
grals involving polynomials of Bessel functions. When 
kp ~ the RKKY interaction decays monotonically as 
l/i?3. When /cpi? ^ 1 all the coefficients have oscil- 
latory behavior and the leading terms decay as l/R^; 
this behavior is characteristic of an ordinary 2DES^ii^, 
except for the fact that B,C,D depend on kp- Density- 
dependent amplitudes are ubiquitous in graphenoi^ as 
well, where nonetheless C ^ D = 0. Finally, we ver- 
ify that {A + C)/A2 = B/X\ for any kp, which is also 

satisfied in a 2DES with Rashba SO couplingi^. 
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FIG. 1: (a) RKKY coefficients between two impurity spins in 
the kpR > 1 regime (Eq. {TSj). We take fc|./(167r^u) = 1. At 
large distances A is parametrically smaller than B,C,D. (b) 
Two ultrathin ferromagnetic strips deposited on the surface 
of a 3D TI (viewed from the top), (c) The integration vari- 
ables in Eq. psp . (d) Interlayer exchange coefficients lij for 
two uniform ferromagnetic monolayers, as a function of the 



interlayer separation. 
kpl = 1000. 



We took 16tv vslSu = 1, kpw = 25, 



We now consider two uniformly magnetizeS^ ferro- 
magnetic monolayers depicted in Fig. [T]3,c. The inter- 
layer exchange coupling per unit length is 



i,j£{x,V,z} 
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where ^l{r,) is the direction of magnetization in the left 
(right) ferromagnet. 
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and Si(fl) is the area per spin in the left (right) magnet. 
When kp = Q and I ^ w ^ d we find a ferromagnetic 
interlayer coupling with strong in-plane anisotropy: 



^xc(/J — 0) — Ixx^'lS^^ + Izz^l^r'^ 
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For kpR ^ 1 we evaluate 7^ numerically and arrive at 



(17) 

with the coefficients displayed in Fig. [TJi. Once again 
there is a strong in-plane exchange anisotropy (/yy ~ 0). 
The last term of Eq. (|17p is a Dzyaloshinskii-Moriya 
(DM) interaction^^ that favors non-coUinear magnetiza- 
tion configurations. The DM vector D = Izxi) alternates 
sign as a function of the separation between the ferro- 
magnets and may be tuned by a gate voltage. 

RKKY interaction in presence of an in-plane electric 
field — It is well-known that the amplitude of inter- 
layer exchange coupling between two magnets changes 
when the non-magnetic spacer is driven out of electric 
equilibrium^^. It is less known that electric fields can 
also induce unconventional, non-Heisenberg types of ex- 
change coupling, provided that the spacer contains spin- 
orbit interactions. Due to its strong SO interaction, the 
surface of a 3D TI constitutes an ideal spacer where such 
an effect might be observable. 

We compute the change in the static spin-spin suscep- 
tibility under a chemical potential yu(r) which is constant 
in time and varies slowly in space (Fig. [2]) : 



5x^M^ 0; Q, 0) = Q ■ n,jM(Q, 0) + 0{Q^), 
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where Q is the wave vector of the chemical potential and 
n,-, is given by the Q — limit of 



/ Tr[a*G(k+,c.)aJG(k_,L.)lG(k,w)] 
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k+ = k + q, k_ = k — Q and the factor of 2 stems from 
summing the second diagram in Fig. [2] For simplicity we 
assume that the unperturbed system has the chemical 
potential at the Dirac point; in such case it is possible to 
evaluate Eq. ([T^ analytically. We obtain 



^Xij ^ 5Xxz{kx5jz+5iz6jx)+5Xyz{ky5jz+6iz5jy), (20) 

where 5ij is the Kronecker delta. 
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FIG. 2: Influence of an in-plane electric field on the RKKY 
interaction: Feynman diagrams for Sxij- We pick the contri- 
bution linear in Q, which is associated with the response to 
a uniform in-plane electric field. 



and we have used Q = Q {cos (j)Q, sin (pq). Fourier trans- 
forming from q to R , we see that the in-plane electric 
field produces non-collinear RKKY interactions between 
two magnetic impurities: 



(22) 



(R X eE), 
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where SJij = —Xfjdxiji'R) and E = — iQ/i/e is the 
electric field. SE^kky = when E||R. The ampli- 
tude of 6EjiKKY equals that of E'rkky (Eq. (HH)) when 
E _L R and £" ~ = hv/{eR^). For i; 5 x lO^m/s, 
i?c [mV/nm] ~ 330/i?[nm]^. Thus moderate electric 
fields suffice to modify the equilibrium RKKY interac- 
tion substantially, especially when the impurities are sep- 
arated by a distance > 5nm. 

For the system shown in Fig. [T|d, the change in the 
interlayer exchange coupling due to an in-plane electric 
field is 
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where we have assumed Z 3> w 3> d. SS-^c ~ when 
rif^ = rifj = and (5£xc may be nonzero even when 
Ex = ^ Ey, i.e. when there is no potential difference 
between the ferromagnets. Moreover, the amplitude of 
S8^c becomes comparable to that of S^c (Eq. (HH)) when 
E > hv/eiD^. The task of evaluating the kp ^ Q coun- 
terparts of Eqs. and is cumbersome and will be 
left for future work. 

Coupled spin-charge dynamics — Eq. ^ implies that 
a time-dependent Zeeman field U will induce a charge 
density n\-a_d on a T- invariant helical metal. For kp ~ Q 
we find 
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which is valid for an arbitrary Zeeman field. Motivated 
by the possibility of experiments with single spins de- 
posited on the surface of a 3D TI we concentrate on a spa- 
tially localized perturbation, whereby t/(q, Vt) = U{H) 
and 
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It follows that 
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The induced current defined via V • ji 



'dtriind is then tardation effects") and assuming kpr ^ 1 we arrive at 
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where we have used the boundary condition jind(?' ^ 
oo) = 0. The charge and current densities at time t are 
determined by the external perturbation at a retarded 
time t — r/v, which reflects the relativistic nature of Dirac 
fermions (with v as the "speed of light" ) . 

For concreteness, consider a single impurity whose spin 
is pointing along x in equilibrium. Under a constant mag- 
netic field B = Boy, the impurity magnetic moment pre- 
cesses and U(t) = X^^[cos{jBot)x + sm{'yBot)z], where 7 
is the gyromagnetic ratio and we have neglected magnetic 
damping. Using 7 ~ 1.76 x 10"Hz/T, v = 5 x lO^m/s 
audi Am ~ 50meVnm2, Eqs. ^ and ^ read 
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where t* = t — r/v and 0{t) = jBot. 
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is large enough 



to be measurable by a field effect transistor located in the 
vicinity (e.g. r ~ lOnm) of the precessing local moment. 

Likewise may be experimentally accessible through 
its associated magnetic field (cx 1/r for large r), to be 
distinguished from the dipolar field (oc 1/r^) originating 
from the impurity moment. The AC characteristics of 
jind distinguish this effect qualitatively from the charge 
pumping phenomenon that occurs in more ordinary fer- 
romagnet / paramagnet interfacesSS. 

Next, we generafize Eqs. and (gHl) to kp ^ 0. 

Linearizing Eq. © in i7 (which is akin to neglecting "re- 



in Eq. ([50]) we have neglected terms oc 
cos(2fci?r), sin(2A:i?r) because their amplitudes are 
small. Hence, when kpr 3> 1, ni^id and jind decay more 
slowly with the distance from the local moment than 
when kp = 0. 

Conclusions. — Our analytical expressions for the cou- 
pled spin-charge response functions in a ballistic helical 
■ . y^mctal underlie a host of novel and potentially observable 
phenomena that arise when magnetic impurities or pat- 
terned magnetic films are deposited on the surface of a 
(29)rpj^ Applied gate voltages and transport currents lead to 
unconventional RKKY interactions between two such im- 
purities or magnetic monolayers, and oscillatory charge 
currents arc generated by magnetic moments precessing 
under an external magnetic field. Similar effects should 
also occur in ordinary 2DES with Rashba SO interaction, 
though in a much weaker form because the Rashba SO 
strength is typically much smaller than the Fermi energy. 
It would be interesting to replicate our calculation for a 
diffusive helical metal, comparing it with Ref. 
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